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Generic suppression of conductance quantization of interacting electrons in graphene 

nanoribbons in a perpendicular magnetic field 
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The effects of electron interaction on the magnetoconductance of graphene nanoribbons (GNRs) 
are studied within the Hartree approximation. We find that a perpendicular magnetic field leads to a 
suppression instead of an expected improvement of the quantization. This suppression is traced back 
to interaction-induced modifications of the band structure leading to the formation of compressible 
strips in the middle of GNRs. It is also shown that the hard wall confinement combined with 
electron interaction generates overlaps between forward and backward propagating states, which 
may significantly enhance backscattering in realistic GNRs. The relation to available experiments 
is discussed. 
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Conductance quantization in quantum point contacts 
(QPCs) and quantum wires represents a hallmark of 
mesoscopic physics [l|, l2j . At zero magnetic field this ef- 
fect can be understood within a noninteracting electron 
picture as quantization of the transverse electron motion 
where, according to the Landaucr-Buttikcr formalism, 
each propagating mode contributes with the conductance 
quantum Go = 2e 2 /h to the total conductance [l|, |2j. In 
a perpendicular magnetic field B the propagating states 
acquire qualitatively new features gradually transform- 
ing into edge states as B is increased |2]-l5j]. Since the 
left- and right-propagating edge states get localized in 
transverse direction at opposite wire edges in sufficiently 
strong magnetic fields, the coupling between them can 
be exponentially small. This, in turn, leads to a strongly 
suppressed backscattering and hence to a drastic im- 
provement of the conductance quantization[2]-|6|]. Taking 
electron interaction and screening in high magnetic fields 
into account leads to new features such as formation of 
compressible and incompressible strips [7(, which are es- 
sential for an interpretation of various magnetotransport 
phenomena in conventional QPCs and quantum wires de- 
fined in two-dimensional electron gases (2DEGs) |7J, |8[ . 

The isolation of graphene [9[ has immediately inspired 
the search for conductance quantization in graphene 
nanoribbons (GNRs). However, in all experiments re- 
ported so far conductance quantization at B = is ab- 
sent [lfj or strongly suppressed [ll| , which by now is well 
understood and attributed to the effects of impurity scat- 
tering and/or edge disorder [12J . In analogy with conven- 
tional QPC structures one would thus anticipate a drastic 
improvement of the conductance quantization in GNRs 
in the edge state regime due to the expected suppression 
of backscattering |J|. Surprisingly enough, the magne- 
toconductance measurements on GNRs reported so far 
show no evidence of the expected improvement of the 



conductance quantization I13l . Il4ll . Even relatively large 
graphene strips (> l^jm) [la . Il6l | do not exhibit quanti- 
zation plateaus at high magnetic fields of high quality as 
routinely seen in corresponding conventional hctcrostruc- 
tures [fj. 

In the present paper we study the magnetoconductance 
of GNRs taking electron interaction on the Hartree level 
into account. Contrary to expectations based on the con- 
ventional edge-state picture of noninteracting electrons 
[4J we find that application of a magnetic field leads to 
a suppression instead of expected improvement of the 
conductance quantization. This behavior is related to a 
drastic modification of the GNR band structure by elec- 
tron interaction leading, in particular, to the formation 
of compressible strips in the middle of the ribbon. These 
features are generic in GNRs, but in contrast to most 
of the distinct properties of graphene [18[ they are not 
caused by the Dirac-likc energy dispersion but rather by 
the hard- wall confinement. 

We consider a GNR attached to semi-infinitive leads 
acting as electron reservoirs and subjected to a perpen- 
dicular magnetic field B, see inset to Fig. [T] The 
ribbon of width w = 50 nm resides on top of a Si02 
insulating substrate (ey = 3.9) of thickness d = 300 
nm, below which a metallic gate is located. The sys- 
tem is described by the standard p-orbital tight-binding 
HamiltonianfH [S| 
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where the summation runs over all sites of the graphene 
lattice, A includes the nearest neighbors only, i r , r +A = 
to exp(i27r0 rr+ A/</ ) o) with to = 2.77 cV, (f>o = h/e being 
the magnetic flux quantum and </> r , r +A = / A ■ dl 
with A being the vector potential. We use the Landau 
gauge, A = (—By,0). The interaction among the extra 
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FIG. 1: (Color online) Conductance of the CNR as a function 
of filling factor for interacting and noninteracting electrons at 
temperatures T = 20 K and 50 K in a magnetic field .8=30 
T (corresponding to Ib/w ~ 11). The arrows indicate the 
filling factors for which the corresponding band structures are 
shown in Fig. J5J . Inset: sketch of the sample geometry. An 
armchair CNR of width w — 50 nm is located on top of an 
insulating Si02 layer (e r = 3.9, thickness d — 300 nm) and a 
gate electrode. 



charges of the density n(r) is described within Hartree 
approximation 
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where the first term describes electron interaction within 

the ribbon while the second term takes the presence of 
the metallic gate on the basis of the image charge method 
into account. The band structure, the potential pro- 
file, the charge density distribution are calculated self- 
consistently using the Green's function technique (see 
Rcfs. [20j, [2l| for details). The magnetoconductance 
through the nanoribbon in the linear response regime is 
given by the Landauer formula 



G(E F: B) = ^-Jt(E,B) 
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where Jfd(E — Ep) is the Fermi-Dirac distribution func- 
tion and Ep denotes the Fermi energy. For an ideal 
system (without scattering), the total transmission co- 
efficient T(E, B) is equal to the number of propagating 
states, T(E,B) = N prop , such that the conductance is 
simply proportional to N prop weighted by fep which is 
different from zero in an energy window ~ AirksT. 

Figure [1] shows the conductance of the ideal nanorib- 
bon for a representative magnetic field B = 30 T as a 



function of the filling factor v = (n) 4>q/B for two repre- 
sentative temperatures, with (n) being the electron den- 
sity averaged across the ribbon. Here, v is tuned by vary- 
ing the gate voltage V g which is applied vs. the grounded 
nanoribbon and thus tunes the electron density. The ra- 
tio of GNR width to magnetic length, w/Ib ~ 11, is 
chosen in accordance with typical experiments [13l Il4j . 
It is important to emphasize that the obtained results 
remain practically unchanged when the system is scaled 
by, e.g. increasing w while simultaneously reducing B 
such that the ratio w/Ib remains constant. In order 
to highlight the role of electron interaction, we compare 
our self-consistent calculations with a noninteracting pic- 
ture. The calculated conductance shows a striking differ- 
ence between the interacting and noninteracting cases. 
First of all, at a given filling factor, the conductance of 
the interacting system is always larger than that one of 
the corresponding noninteracting system. Second, the 
perfect quantization steps calculated for the noninteract- 
ing picture are destroyed as the interaction is turned on, 
and the conductance develops pronounced bump- like fea- 
tures. Note that the elevated temperature smears the 
conductance bumps to some extent, but they still dom- 
inate the conductance even at T = 50/f. We also note 
that we performed the magnctotransport calculations for 
a high-k material (e r = 47) and a gate closcby, d = 5 nm 
when the electron interaction is strongly screened (not 
shown here) . We find that even in this case the bumps are 
suppressed but still clearly dominate the conductance. 

We proceed by interpreting the suppression of conduc- 
tance plateaux in terms of interaction-induced modifica- 
tions of the energy dispersion. The evolution of the band 
structure as a function of v in the interval covering a 
representative bump, 2.9 < v < 5.7 (corresponding to 
arrows (a)-(c) in Fig. [T|) is presented in Fig. [5] both for 
interacting and non-interacting cases. The dispersion re- 
lation for non-interacting electrons shows flat regions cor- 
responding to the Landau levels in bulk graphene[18ll23| 
and dispersiveness states close to the GNR's boundaries 
representing familiar edge states [5j. Note that the po- 
sition and the number of propagating states at a given 
energy are determined by the intersection of the Fermi 
energy level with the corresponding subbands. 

For non-interacting electrons changing the gate voltage 
results in a shift of the Fermi level but does not modify 
the subband structure. Qualitatively new features arise 
when the- electron interaction is taken into account. One 
of the most distinct features is that the dispersionlcss 
state in the center of the GNR (corresponding to the 1st 
Landau level (LL)) gets pinned to the Fermi energy thus 
forming a compressible strip. These strips are marked in 
yellow in Fig. |2](b)-(d); following Suzuki and Ando J24| 
we define a compressible strip as a region where the dis- 
persion lies within the energy window \E—Ep\ < 2irkBT. 
The compressible strips form because in the above energy 
window the states are partially filled (i.e. < fpD < 1) 
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FIG. 2: (Color online) Evolution of the band structure of the GNR at different filling factors corresponding to arrows (a) - 
(d) in Fig. [l] Left and right parts of the panels correspond to the interacting and non-interacting case, respectively. In order 
to align noninteracting and Hartree bands the one-electron dispersions have been shifted along the energy axis by the average 
Hartree energy. Gray fields mark the energy window [Ef — 2irkbT,EF + 2nkbT]; yellow fields mark the compressible strips. 
The dotted line shows Ef- The black full circles mark the intersections of the Fermi level with the dispersion curves, thereby 
identifying the propagating states at Ef- In (a) the dispersionless states are marked according to the corresponding Landau 
levels (LL) of the bulk graphene. 



and hence the system has a metallic character. Due to 
the metallic behavior, the electron density can easily be 
redistributed in order to effectively screen the external 
potential [7| . The compressible strips can form only if the 
confining potential is sufficiently smooth [7[. The GNRs 
have a hard- wall confinement and hence the compressible 
strips can form only in the center but not for the edge 
states. Note that the existence of compressible strips in 
graphene has been recently demonstrated by Silvestrov 
et oZ.Q. 

Because of the pinning of the LL to the Fermi energy, 
changing of the filling factor leads to a significant distor- 
tion of the dispersion curves. For a given B, the larger 
the gate voltage (and therefore v) is, the stronger the 
bands arc distorted in comparison to the noninteracting 
picture (c.f. (a)-(d) in Fig. [2] ). This distortion even- 
tually leads to the bumps in the conductance. Indeed, 
according to Eq. ^ the conductance is given by the 
number of propagating states averaged in the energy win- 
dow \E — Ep\ < 2&bT '. For noninteracting electrons the 
dispersion relation is not changed as v is varied and the 
number of propagating states remains always the same, 
Np r0 p = 3 (see right panels in Fig. [5] (a)-(c)). This, 
according to Eq. ([3]), leads to a conductance plateau 
G = 3Go- In contrast, for interacting electrons the dis- 
persion relation gets distorted and there is always an en- 
ergy interval in the window \E — Ep\ < 2ksT where the 
number of propagating states exceeds that for the non- 
interacting case. This is illustrated in the left panels in 
Fig. E](b)-(c) for E = E F where N prop — 5. As a result, 
the conductance exceeds its noninteracting value of 3Go 
exhibiting the pronounced bumps. 

With further increase of v, the compressible strips 
pinned to the Fermi level form not only in the center 



of the strip, but further away from the center (as illus- 
trated in Fig. E](d)). Note that the second compressible 
strip in Fig. [5] (d) leads to the formation of a bump in 
the conductance in the region 6 < v < 9. 

Let us now discuss in detail a structure of propagat- 
ing states of the interacting electrons in GNRs. Figures 
[3] (a), (b) show the electron density and the confining 
potential for a representative filling factor v = 4.7 [(b)- 
arrow on Fig. ((T|)]. The distribution of charge density is 
highly nonuniform showing charge accumulation at the 
boundaries |2l|,[22|. There are two types of states, which 
have a different microscopic character. The first type 
[1,2,3 states in (a)-(c)] corresponds to edge states propa- 
gating near the boundaries and have the same structure 
for interacting and noninteracting cases. The second type 
[states 4 and 5] corresponds to the states which form com- 
pressible strips in the center of the ribbon as discussed 
above. The most prominent feature of these states is that 
their direction of propagation is opposite to that one of 
the edge states residing in the same half of the GNR. This 
is in contract to the noninteracting picture, where due to 
the presence of a magnetic field, forward and backward 
propagating states are localized at different boundaries 
by Lorcntz forces. This unusual behavior can be inter- 
preted in terms of a semiclassical analogue. The electrons 
scattered at the boundaries are described by skipping or- 
bits. Besides the hard-wall potential walls provided by 
nanoribbon's edges, there are two additional walls origi- 
nating from the self-consistent potential which, together 
with the outer walls of the GNR, form triangular quan- 
tum wells at the ribbon's edges [Fig. E2(b)]. Electrons 
which strike the left side of the right triangular quantum 
well propagate in the same direction as the electrons that 
strike the left edge of the nanoribbon as schematically il- 
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FIG. 3: (Color online) (a) Electron concentration n(y), (b) 
self-consistent potential Vh across the GNR, and (c) the band 
structure at v = 4.7. (d) The square modulus of the wave 
functions at Ef of forward (F) and backward (B) propagating 
states (solid and dashed lines correspondingly) marked in (c) . 
For the sake of clarity the electron densities, potential and 
the wave functions are averaged over two adjacent slices and 
3 adjacent sites of the same slice). Inset in (b) illustrates 
classical skipping orbits. 



lustratedin Fig. [31(b). 

This feature of propagating states in high magnetic 
field makes GNRs much more sensitive to the effect of the 
disorder in comparison to conventional split-gate struc- 
tures defined in 2DEG. Indeed, for interacting electrons 
in GNRs the overlap between the backward (4B, 5B) and 
forward (1F-3F) propagating states is significant. In re- 
alistic GNRs with disorder this would result in a strong 
enhancement of backscattcring, which, in turn, can lead 
to a further distortion of the conductance (in addition 
to bumps that arc present even in ideal GNRs without 
disorder). 

Note that the features of the band structure and char- 
acter of propagating states in GNRs discussed above are 
not caused by the Dirac-likc energy dispersion but rather 
by the hard-wall confinement at the boundaries. These 
features of the GNRs resemble corresponding features of 
cleaved-edge overgrown (CEO) quantum wires [171 ] that 
also have a hard- wall confinement. We therefore expect 
that magnetoconductance of CEO also should exhibit 
suppressed quantization of high field (Note that we were 
not able to find any reports on magnetoconductance mea- 
surements in CEO wires at high magnetic field). 

We continue by relating our results to the available 
experimental data. We are not aware of any studies re- 
porting a drastic improvement of the conductance quan- 
tization in GNRs by perpendicular magnetic field. The 
observed conductance in narrow GNRs exhibit irregular 
14j or bump- like features [13(, and the wider structures 
show pronounced bumps superimposed on conductance 
plateaus [15j, [iff . Even though this is consistent with our 
findings, this can hardly be regarded as a definite exper- 
imental validation of our predictions. We thus hope that 



our work will motivate systematic studies of the magne- 
toconductance that will shed new light on properties of 
interacting electrons in confined graphene systems. 

In conclusion, we have shown that applying a perpen- 
dicular magnetic field to a GNR containing an interac- 
tion electron gas leads to a suppression instead of ex- 
pected improvement of conductance quantization. This 
surprising behavior is related to the modification of the 
band structure of the GNR due to the electron interac- 
tion leading, in particular, to the formation of compress- 
ible strips in the middle of the ribbon and existence of 
counter-propagating states in the same half of the GNR. 
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Swedish Research Council (VR). 
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